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Abstract 

We study a system of two tunnel-coupled quantum dots, with the first dot containing interacting 
electrons (described by the Universal Hamiltonian) not subject to spin-orbit coupling, whereas the 
second contains non-interacting electrons subject to spin-orbit coupling. We focus on describing the 
behavior of the system near the Stoner transition. Close to the critical point quantum fluctuations 
become important and the system enters a quantum critical regime. The large- approximation 
allows us to calculate physical quantitites reliably even in this strongly fluctuating regime. In 
particular, we find a scaling function to describe the crossover of the quasiparticle decay rate 
between the renormalized Fermi liquid regime and the quantum critical regime. 
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I. INTRODUCTION 



The transport of electrons through mesoscopic systems at low temperatures is a coherent 
process. The manifestations of coherent electronic motion are weak localization, Aharonov- 
Bohm oscillations, persistent current, etc.'^ 

Though many mesoscopic effects can be explained in the framework of non-interacting 
electrons, there is a vast amount of evidence showing that collective effects of the electron 
spin are important in predicting the behavior of the system. Apart from its fundamen- 
tal interest, studying interaction- related effects on the electronic spin is important from a 
technological point of view. 

The spin of the electron couples to the external magnetic field and to the orbital degrees of 
freedom. This spin-orbit coupling (SO) is caused by a non-zero electric field in the laboratory 
reference frame that is transformed into a magnetic field in the electron's rest frame. In bulk 
systems the SO coupling results from the absence of inversion symmetry in the crystalline 
latticd^ (Dresselhaus term). In finite size systems, such as metallic grains or semiconductor 
quantum dots, an additional contribution to SO coupling comes from the structure inversion 
asymmetrj^ (Rashba term), the simplest example of which is a two-dimensional electron gas 
(2DEG) confined to an interface, in which the confining electric field perpendicular to the 
2DEG is the source of the SO coupling. 

In diffusive and ballistic/chaotic mesoscopic systems the kinetic term in the full Hamilto- 
nian is well described by Random Matrix TheorjJ^El (RMT). RMT has been very successful in 
describing the ensemble averages of one-particle spectral correlations as well as correlations 
of eigenfunctions. 

RMT describes the universal zero-dimensional limit in a mesoscopic system. Its regime 
of validity is when all time scales (the spin-orbit relaxation time tso and the inverse mean 
level spacing are much larger than ergodic time Terg = h/Ex- Alternatively, all relevant 
energy scales should be smaller that Thouless energy Et (for a diffusive dot of linear scale 
L, Et — hD/L^, where D is the diffusion constant, while for a ballistic/chaotic dot Et — 
hvp/L). 

Even though tso ^ Terg defines the universal limit, to decide if the SO coupling is 
important for a particular physical process, tso should be compared to other characteristic 
time scales.'^'^ As the SO coupling is increased from zero in a noninteracting system, its 
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effects begin to become important for physical quantities when the inverse SO relaxation time 
is comparable to the mean level spacing 6 ~ t^q. In interacting systems near a degeneracy 
point between two ground states of different spin, even much tinier SO couplings can have 
order one effectsP 

RMT systems can be classified according to the presence or absence of time reversal 
(TR) and spin rotation symmetries. They fall into three major categories described by 
the "classical" RMT ensembles introduced by Dyson. The systems with both symmetries 
preserved belong to the Gaussian Orthogonal ensemble. Systems with broken TR symmetry 
(e.g. by an external magnetic field) are described by the Gaussian Unitary ensemble. Finally, 
systems with TR preserved and broken spin rotation symmetry belong to the Gaussian 
Symplectic ensemble. 

In the universal regime, the SO coupling has significant effect on spectral properties of 

eigenvalues and eigenfunctions. One can relate the spin-orbit scattering length L^o to a SO 

crossover energy scale Ex = Ej-(^j^Y ■ energies below Ex, the one-particle term in the 

Universal Hamiltonian is modeled by a Gaussian symplectic random matrix. If one wants 

to describe energies both above and below Ex, one has to use the RMT ensemble which is 

in a crossover between the GOE and the GSE classes.'^'^ 

The interactions in mesoscopic systems at low temperatures are described by the Universal 
HamiltonianPinillil 

Hu = Yl e-cj,,,c„,, + ^N' - JS^ + XT% (1) 

a,s 

where N is the total particle number, S is the total spin, and T = ^C/3,j^C/3j. The 
Universal Hamiltonian contains a charging energy U (direct channel), a Stoner exchange 
energy J (spin channel) and a reduced superconducting interaction A (Cooper channel). A 
renormalization group (RG) analysis reveal^^^*^ that this is the low-energy effective theory 
for weak coupling, although other effective theories and other ground states can be accessed 
for strong coupling.'^^ 

For small normal metallic grains and non-superconducting quantum dots with a fixed 
number of particles the exchange interaction is the main contribution to electron-electron 
interactions. The short range part of electron-electron interactions causes the ferromagnetic 
Stoner instability at large values of exchange energy J. 
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In the absence of SO coupling the total spin of the system and its z-projection 5*^ 
commute with kinetic energy term and are good quantum numbers. Typically for metallic 
grains the exchange constant J < 5. 

For weak exchange interaction J <^ 5 the spin of the ground state for odd number of 
electrons is 1/2. As J gets larger, there is a non-zero probability to obtain a ground state 
with S > 1/2. This happens when the cost in orbital energy to promote an electron to the 
next level is less than the energy gain due to the exchange interaction. As J approaches (5, 
the total spin of the system grow^, and at J = 5 the system undergoes a phase transition 
(the Stoner transition) into a "bulk" ferromagnetic state. For J > 6 the magnetization of 
the system is proportional to the number of electrons N. 

In the presence of SO coupling the total spin does not commute with full Hamiltonian 
[H, S^] 7^ 0. While the dominant effect of the electron-electron interaction is to organize the 
states according to total spin S, the SO term produces matrix elements between states of 
different spin, which randomizes spin, and also leads to sample-to-sample fluctuations of the 
matrix elements of the electron-electron interactioiJ^ and the suppression of the exchange 
interaction.'^'^ 

When J ^ 750, the SO coupling is simply ignored (unless one is near a degeneracy 
between ground states of different spin.l^ In the opposite limit J <ti 6 and J <ti "jso electron- 
electron interactions are suppressed and expectation value of total spin in ground state 
(S) < 1/2. The interesting regime is when 750 ^ J ^ ^- In this case the exchange 
interaction is not completely suppressed, and the fluctuations of total spin are comparable 
to its expectation value. This regime is driven by the combined effect of spin-orbit scattering 
and electron-electron interactions. 

We study the regime where the system is near the Stoner instability J ^ 6~ . If SO 
coupling is absent, there are no quantum fluctuations of the spin, and one obtains a se- 
quence of metamagnetic transitions with the true Stoner transition being the accumulation 
point.'^In the presence of SO coupling, at low energies, the behavior of the system is domi- 
nated by quantum critical fluctuations leading to the formation of a quantum critical regime 
(QCR).™ 

Imagine that one is at some J < 6, but that 1 — J/5 ^ 1. Even close to the transition one 
can think of two different regimes of energy separated by a many-body crossover scale Eqcx, 
which will turn out to be simply related to the single-particle RMT crossover scale Ex^For 
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uj <^ Eqcx the system behaves as though it were a renormahzed Fermi hquid, with altered 
Fermi hquid parameters and a quasiparticle decay rate going as u;^. ISM] ^j-^g other hand, 
for Eqcx ^ ^, the behavior is controlled by the quantum critical point. The change of 
behavior as one increases u is described by a universal scaling function F{uj / Eqcx)- 

The critical point and QCR are dominated by many-body quantum fluctuations, and 
thus the scaling functions cannot be calculated perturbatively. However, it turns out that 
as long as Eqcx,^ ^ one can use a large-A^ approximation with min(^^^^, j) playing 
the role of the large A^.^^ This allows us to compute the scaling functions reliably. 

From the point of view of experiment, the key point is that one can control Eqcx, which 
is a many-body scale, by tuning a single-particle crossover energy scale Ex- Thus, at a 
fixed value of the parameter J, one can tune oneself into and out of the QCR by tuning a 
single-particle knob. 

As a prerequisite to describing the system near Stoner transition, we consider the non- 
interacting case and calculate ensemble-averaged one and two particle Green's functions 
for electrons in the first dot coupled to the second dot in crossover between GOE and GSE 
ensembles. The one particle Green's function is unchanged by crossover, though it is modified 
by interdot coupling. The two particle Green's function is the sum of the contributions due 
to diffusorp2lmode and Cooperon modes. Bot contributions depend on the ratios of crossover 
parameter Ex2, interdot coupling parameter Eu, and measurement energy uj. 

It may seem counterintuitive that one can use non-interacting wavefunction averages to 
describe the behavior of a system with strong many-body fiuctuationd^^IIHIISEl^ ^^^^ ^j^jg 
goes hand in hand with the use of the large- approximation. This is because there is no 
wavefunction renormalization to leading order in the the large- A^ approximation. 

II. MODEL DEFINITION 

We consider a system of two quantum dots (metallic grains) coupled to each other by 
tunneling (see Fig. [T]). The motion of electrons can be either diffusive or ballistic/chaotic: 
in either case the single-particle energies and wavefunctions are controlled by RMT, which 
is all that we require. 

For the non-interacting system the ensemble-averaged spectral and eigenvector corre- 
lations can be computed by RMT."^ The first dot belongs to the GOE, since it has no 
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Dot 1 

with Stoner interaction 
no spin-orbit coupling 




Dot 2 

with spin-orbit coupling, 
no Stoner interaction 



FIG. 1: Two coupled quantum dots in GOE-to-GSE crossover. 

InSb 




GaAs 



Insulating layer 



FIG. 2: (Color online) The vertical arrangement of the dots allows us to get rid of the charging 
energy. 

spin-orbit coupling. The second dot has weak spin orbit coupling that drives it into the 
GOE— s-GSE crossover, characterized by a crossover scale Ex2- We also assume that there 
is a Stoner exchange interaction in the first dot. No interactions are present in the second 
dot. The tunneling between the dots gives rise to another crossover scale Eu, where Eu/S 
is the dimensionless conductace between the two dots. 

In Fig. |2] one can see a more realistic picture of the system. In an experimental setup 
the lower dot could be made of GaAs (with significant exchange interaction but tiny SO 
coupling), while the upper dot could be made of InSb (with large spin-orbit coupling). The 
choice of the vertically coupled geometry will be discussed below. 

In the low-energy limit interactions are described by Universal Hamiltonian.^^ ^ * ^^ * ^^ ' ([T]) 
For our system Hu has the form: 
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H=Y1 H^A,sC,o,s 'J^'+Y. HSloio,sC.o,s + Yl y^o,oi4o,sC,o,s + H.c 

iojos 



iojos t^oi^os iofJ-OS 
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where H^"^^ contains the effect of spin orbit couphng in the second dot. In ^ we have omitted 
the superconducting term as irrelevant to our model. We also choose a vertically coupled 
geometry for our system to minimize the change in charging energy when the electron hopes 
from one dot to another-f^ If this energy is smaller than all other relevant scales, then 
the charging term can be omitted (or absorbed into the chemical potential) since the total 
number of electrons in the two-dot system remains unchanged. 

The (/X, r) label the basis of the two coupled quantum dots without interaction, that is, 
it is the set of eigenstates of H^^^ + H^'^^ + V. Here /i is the orbital quantum number and r is 
a twofold degenerate Kramers index. In this basis the a-th component of total spin read^ 



S'' = Y 4j-fc^os' = E iMXAr^.r', (3) 

where io is an Orthogonal basis in the first dot alone (it could be the eigenbasis of H^^^ but 
it does not have to be), and the matrix element is defined as: 

i^Xl' = Yr,M<^^s)^^^A^o,s'). (4) 

iOySS 

We reiterate that the first summation in Eq.([3]) and the summation in Eq.(|4]) is over an 
Orthogonal basis in the first dot alone, while the second summation in Eq.(|3]) is over the 
eigenbasis of the total non-interacting Hamiltonian H^^^ + H^'^^ + V. Also, a" / are the Pauli 
matrices, ipuriio, s) is the wave function of the state /i, r in the first dot. 

We use Eq.([2]) to calculate the partition function Z = Tr[exp(— /3if)], using the imaginary 
time path integral formalism: 

Z = Trie-'''') ^ Z= fH Vc^rVc.rVh e" ^'^^ (5) 
where the Euclidean Lagrangian is: 
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|2 



Ihl 

1^ = -jj + c^r{dt + e^)c^r - h ■ S. (6) 

In Eq.([5]) we used the Hubbard-Stratonovich transformation to decouple the interaction 
at the expense of introducing an additional bosonic field h representing the order parameter. 
The Cfj,r and c^r are Grassmann variables. 

After switching to the Fourier representation the fermionic fields c, c are integrated out. 
The resulting action for h is expanded to second order to obtain: 

^^// ^ E [j - Ex,,Eu)] (7) 



n,a 



6^ 6,y lUJr] 



flT,l/T 

where a;„ = 2nn/P, J = J/(5i is a dimensionless exchange constant, Np{e^) is the Fermi- 
Dirac occupation of the state /i, and 6i is a mean level spacing for the first dot. 

Deep into the crossover Eu, Ex2 > 5i, we replace jM'^p''^ by its RMT ensemble average. 
This is justified because in the limit when ^ — > oo the spectral average on a single 
sample is the same as the ensemble average. The corrections to this vanish in the large- 
limit. This is one of the ways in which we use the large- approximation. 

The relevant four wavefunction correlator hidden in |M°p) is calculated in Appendix [b| 
We also replace the summation over energy eigenstates by energy integrations. Assuming a 
constant density of states we obtain 



fn{(3, Ex2, Eu) 



El - El 



+ Ex^Eu - El ^ E\^ + Ex^Ey - E^ 

El + \Un\ ^2+ k„| J 

where the interdot tunneling energy scale Eu, the SO crossover energy scale Ex2 in the sec- 
ond dot, and the energies £'12 (which are functions of Eu and EX2) are defined in Appendices 

The instability point is obtained by setting fo^P, Ex2, Eu) = J^^ ■ For the coupled-dot 
system the quantum phase transition takes place at J = 1, or J = 5i, the same result as for 
one uncoupled dot independent of the crossover energy scales. 

We investigate the limit when Eu ^ Ex2 ■ In this limit Eu is the only relevant parameter 
that controls both the coupling between dots and the degree to which spin rotation symmetry 
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is spoiled in the first dot. In this hmit the scahng function becomes /„ = Eu/{Eu + \ujn\)- 
Close to the transition the smallness of 1 — J allows us to introduce a new scaling function 
Fn that describes the interacting system near Stoner transition. The effective action now 
becomes: 

^ n,a 



JEu V Eqcx 



The scaling function F„ in (11) describes how a physical quantity behaves when one goes 
from the renormalized Fermi liquid to the quantum critical regime. The new characteristic 
energy scale Eqcx = Eu{l — J) can be used to tune the system into the QCR. By changing 
the single particle parameter Eu in Eqcx one can access the QCR governed by interactions. 

The Fig. [3] shows the phase diagram in {uj,T) vs. J coordinates. For J > Si, the system 
is in a "bulk" Stoner phase where the magnetization is proportional to the volume of the 
system. When J < 61 and the measurement energy u satisfies the inequality Eu <^ < Et 
one enters an approximate spin-rotation-invariant universal regime described by Universal 
Hamiltonian Hu. Here the total spin of the system is (approximately) a good quantum 
number. Lowering the energy to ~ Ejj brings us to regime where the system starts seeing 
spin-orbit coupling and spin fluctuations become important. Below the line E = Eu{l — J) 
is the renormalized Fermi liquid regime. Above this line, away from the Stoner instability 
point, is the non-universal regime that depends on many parameters. Close to the instability 
point J = 5i is a Quantum Critical Regime controlled by a single parameter Eqcx- One 
can change the single-particle parameter Ejj in Eqcx to access the many-body regime. 

We proceed to calculate the quasiparticle decay rate near the critical point. Since the 
particle decays by interacting with quantum fluctuations of collective spin, the decay rate 
can be obtained from the spin-spin correlation function {S"'{t)S''{t')) which, in turn, can be 
measured by NMR or EPR. In Fourier space the S"'S^ correlator can be expressed through 
the bosonic field h"" as follows (see Appendix [C] for more details) 

{S''itUn)S\-tUJ^)) = + ^{h''itUJn)h\-tUJn)). (12) 

Calculating the {h°'h^) correlator 



9 



E,T 

Et 



SR Invariant H,j 
Regime 




Bulk 
Stoner 



^1 



FIG. 3: (Color online) Phase diagram in the E vs. J showing different regimes. 



one obtains the following spin-spin correlation function 



^JEr 



(13) 



Oab 

2J 



1 



2Er. 



(14) 



Switching back to the real time formalism {iun ^ uj + iri,ri O'^) in Eq.(14) one obtains 
the spectral function of spin excitations 



-2^ 



S''{lj)S\-uj) 



6. 



2Ei 



u ^ 



ab 



(15) 



The graph of the spectral function ( 15 ) is shown on Fig. 4 



The decay rate of quasiparticles is found by estimating the lowest-order interacting self- 
energy diagram I ^ with interaction V: 



V = -hS = - V hM^" 



(16) 



flT,VT 



The imaginary part of self-energy T,^^^ is evaluated to 



= 5, 



ab 



J 

16n 



In 



Efie^ + E^^^) 



(17) 
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FIG. 4: (Color online) Spectral function for spin-spin excitations. 
The decay rate F for various regimes is plotted on logarithmic scale in Fig. |5j 

III. CONCLUSION 

In this paper we have studied a system of two tunnel-coupled quantum dots (small nor- 
mal metallic grains) near the Stoner transition of the first dot. The first dot has interacting 
electrons but no spin-orbit coupling, while the other has spin-orbit coupling, but no inter- 
actions. The two single-particle crossover energies are Eu, which measures the tunneling 
strength between the dots, and which measures the spin-orbit crossover scale in the 
second dot. Electrons tunneling between the dots carry information about spin-rotation 
invariance breaking to the first dot, and produce quantum fluctuations of the first dot's 
spin. 

Our focus is on the regime near the Stoner transition when the exchange interaction J 
is comparable to the mean level spacing 6i in the first dot. This regime is characterized by 
quantum critical fluctuations rising from the interplay between the spin-orbit and interaction 
parts of the Hamiltonian. For this Quantum Critical Regime we derived the scaling function 
describing the behavior of system observable near instability point J = 5i as a function of 
the measurement energy u. The scaling function itself is dependent on a single dimensionless 
ratio uj/Eqcx, as opposed to other parameter regimes where a physical property can depend 
separately on all the energy scales Eu, Ex2,^, 

As an illustrative example we compute the scaling form of the quasiparticle decay rate, 
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FIG. 5: (Color online) Quasiparticle decay rate (solid line) on a logarithmic scale for different 
regimes. The dashed line represents In uj. Panel (a) shows the decay rate in quantum chritical 
regime. The decay rate in the regime of renormalized Fermi liquid is in panel (b). Panels (c) and 
(d) show the decay rate in non-universal and Universal regimes respectively. 

which can be measured by nonlinear conductance measurements. It has a Fermi liquid-like 
form for uj <^ Eqcx, with the decay rate going as ( eqcx )'^' However, for u ^ Eqcx it goes 
as log{uj /Eqcx)- 

One of the main conceptual points we wish to make is that there is an intimate relation 
between the single-particle crossover energies and the many-body quantum critical crossover 
scale. In the simplest case Eu -C Ex2 this relation is Eqcx = Eu{l — J/S). Access to the 
quantum critical regime can be tuned by changing a single-particle parameter. 

An important open question is the effect of quantum criticality on Coumlomb Blockade, 
that is, how are the distributions of the peak positions, heights, and widths affected by 
quantum criticality. We hope to explore this and other issues in future work. 
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APPENDIX A: GOE TO GSE CROSSOVER IN THE SYSTEM OF TWO COU- 
PLED DOTS 

In this appendix we derive one and two particle Green's functions for two coupled dots in 
crossover between GOE and GSE ensembles. To reduce complexity we consider less general 
(but relevant to our system) situation when only second dot is in GOE to GSE crossover. 
The first dot and the hopping bridge belong to GOE ensemble. The generalization where 
all parts of the system are in crossover can be worked out without difficulty. 

The derivation will be rather sketchy for the full derivation in case of GOE to GUE 
crossover can be found in Ref. p]. 

The Hamiltonian (kinetic part) of two coupled dots is: 

(M) 

where H12 are the Hamiltonians for dot 1 and 2, V describes coupling between two dots. 
Following RMT one considers the elements of H12 and V as Gaussian random variables 
(quaternions) with zero mean. In the crossover between GOE and GSE Hamiltonians H12 
take the form: 



H, = 



(A2) 



Similarly, 



V = 



V^^ (g) + Vy^ (^Ty + (g) T, 



, (A3) 

VI + 3r2 ^ ' 

where ifg and H'^.y^ are real symmetric and real antisymmetric matrices. and V^y^ are 
real and imaginary parts of quantum matrix V (note that elements of and V^y^ are real 
numbers). The Tj matrices are related to Pauli matrices as = iufc, k = x,y, z. 
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The Xi and F are crossover parameters. The denominators in Eqs. (A2) and (A3) keep 
mean level spacing constant when Xj and F change. 

In calculations below we assume Xi = F = 0, so the first dot and the bridge belong to 
GOE; the mean level spacing 5i = 82 Ni = N2. 

The elements of Hi and H2 are independent random variables with correlations between 
symmetric and antisymmetric parts 

) = ^ ( ± ^-^^n. + s^,5^,) (A4) 

where indices s(a) stand for symmetric (antisymmetric); Ni is the size of matrix Hi and 5i 
is the mean level spacing (we assume that Ni = N2, which means 61 = 62). Correlation 
between full matrix elements in crossover is 



1 + 3X2 



1 + 3X2 



(A5) 



Here C,i is a "spin" index that numerates elements of r matrices. 
For V matrix correlations between matrix elements are 



NidW 



(A6) 



where primed (unprimed) indices belong to the second(first) dot; U is a dimensionless pa- 
rameter controlling coupling between dots. 

One particle Green's function for coupled dots is: 



G= {E®I-H)'^ 







) H 


' Gil G12 1 




I G21 G22 / 



(A7) 



E-Hi -V 
-V^ E-H2 

Following the steps in Ref . one can obtain the system of Dyson equations for RMT av- 
eraged Green's functions Gu and G22- In large-N approximation only the rainbow diagrams 
contribute. In the limit of weak coupling the solution for Gu is 



{Gab,!) ^ — ^ab^iiaii 



NlSi 



(A8) 
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where dimensionless energy e = 2^^. 

The two particle Green's function in the first dot can be found from the system of Bethe- 
Salpeter equations.'^ This system describes contribution of ladder and absolutely crossed 
diagrams. 

Expression for the full two particle Green's function is 

{Gf^^,{E)Gt,{E + u)) =D, + C, (A9) 
Contribution of ladder diagrams Di is 



Dl = 5ad5i 



be 



27r 1 

Nf6i -iu 



l + i 



Ex2+Eu 



1 Ex^El 

2 0^3 



Contribution of absolutely crossed diagrams is 



(AlO) 



bd 



27r 1 r n++nri , , ^ n+ -nri 



where 



n 



1 + ,:^^^ 



nri 



1 + 2^)2 + ^ 



(All) 



Crossover energy scales Ex^ and Eu are defined as Ex^ = SX2N161/7!' and Eu 
2UNi6i/7i. 



APPENDIX B: CORRELATION OF FOUR WAVE FUNCTIONS 

Consider the matrix element average (|M"p). More generally, 
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where (i, s) is the basis of first uncoupled dot, and (m, r) is the basis of two coupled dots 
without interaction; ipmriis) is the wave function of electron in the first dot in (m, r) basis. 

In full analogy with derivations in Ref. for GOE to GSE crossover one gets the following 
expression for four wave function correlator 
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(B2) 



where 3?[-Di] and 3?[Ci] are real parts of diffuson and Cooperon contributions to the two 
particle Green's function. Mean level spacing 5 of coupled dot system is 5 = 5i/2. 



Comparing indices in Eqs. (Bl) and (B2) after summation over i and i' it is easy to 
see that the diffuson contribution is times larger than that of Cooperon. Therefore, the 
Cooperon contribution is ignored in large-N approximation. 

Substituting (B2) into (Bl) and using results of Appendix [A| for 3?[-Di] one obtains 



\ mr.m r ) 



Si 



u 



rn'r' ,mT) ^"■'^ QO^ Z?2 _ z?2 



327r Ei - Ef 



E\ + Ex^Eu - El E\ + Ex^Eu - E. 



(B3) 



u;2 + El + El 

Here Ejj and Ex2 are the crossover energy scales defined in Appendix |Xj Energy scales £"1^2 
are equal to Ei^2 = cli,2Ex2: where 



«l2 



+ 26 + 1 ± v/(462 + 2b+ 1)2 - 462 



with b defined as 6 = Eu / Ex^ ■ 



APPENDIX C: SPIN-SPIN CORRELATOR 

The spin-spin correlator (spin Green's function) is defined by: 

{S\t)S\t)) = Z-^ J VhVf]Vri S''it)S\t)e-^ (CI) 

where S"'{t) is a component of total spin of the system. Let's split the action S (defined by 
Eq. (|6])) in two parts S = Si + S2. Here 5*1 is the part of the action containing spin S"", 
Si = — J dth"' S"' (t) , and S2 = S — Si contains everything else. 
Then spin-spin correlator can be written as 
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{S\t)S\t)) = j VhVrjVr^ 



d d 



(t) dh\t' 



-e Me 



(C2) 



Integration by parts in Eq. (C2) transfers functional derivative on exp(— 6*2) term. Per- 



forming differentiation one obtains relation 



{S''{t)S\t)) = ^ 5at + 7^ mt)h\t)) 



2J 



4J2 



In Fourier space relation (C3) reads 



(C3) 



1 



(C4) 



APPENDIX D: QUASIPARTICLE DECAY RATE 



The interacting self-energy of electron * 4 in Matsubara formalism is evaluated to 



^2 . . - 77, (Dl) 

[2,71 Y tpn — lUJn — UJ tUJn — ^ 

where Qo and Pq are non-interacting Green's functions for electron and bosonic excitation 
and A{uj') and H{uo") are their spectral representations; ujn = 2n'K / (3 and p„ = (2n+ l)vr//3 
are even and odd Matsubara frequencies. 

After summation over uj^ imaginary part of self energy reads 



c>fv(l)l 



53 [S 



^J\WLj\')H{u:-J)dJ. 



(D2) 



Here we employed the non-interacting expression for electron spectral function A^{uj') = 
ATI assuming small broadening of elergy levels. Spectral function H{uj) for bosonic 

excitations is 
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Using Eqs. (B3) and (D3) we finally obtain the decay rate near the pole 
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